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A new variational perturbation theory is developed based on the q— deformed oscillator. It is 
shown that the new variational perturbation method provides 200 and 10 times better accuracy 
for the ground state energy of anharmonic oscillator than the Gaussian and the post Gaussian 
approximation, respectively, for weak coupling. 
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There have been many proposals to establish methods of extracting nonperturbative information from quantum 
mechanical and quantum field theoretic systems, such as 1/N expansion and methods based on Dyson-Schwinger 
equation . One of the successful ways to obtain nonperturbative information is the variational perturbation theory 0, 
Many variants of the variational perturbation theory based on the Gaussian approximation have been used to 
study various aspects of physical systems 0, 0, IE S IS S ED • Although the variational perturbation theory provides 
one with systematic correction terms to the variationally determined approximation, it has a limitation that only the 
Gaussian wave function(al) may be used as a variational trial wave function(al) for the most physical systems. In 
this paper we attempt to establish a new variational perturbation theory based on the q— deformed oscillator |lllll2| . 
which provides a better approximation than that based on the Gaussian approximation. 

Quantum anharmonic oscillator has been frequently used in developing various approximation methods in quantum 
mechanics and quantum field theory | ll 1 1 Ifi 1(1 1 1 :i| . To develop a new variational perturbation method we will consider 
the anharmonic oscillator given by the Hamiltonian, 

We first illustrate some essential procedures of the variational perturbation theory based on the Gaussian approxi- 
mation. We start by writing the Hamiltonian Q as the sum of the Gaussian Hamiltonian, Hq, and the perturbation 
^3 • term, Vi, as 



H ' where 
C3 



H = H G + Vi, (2) 



H G ^-f + ^\ Vl = ^^^ + h\ (3) 



Since the Gaussian Hamiltonian is that of a simple harmonic oscillator, it can be written in terms of the creation and 



the annihilation operator, Aq = \ — — x H — p, in the quadratic form: 
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Moreover, the Lie algebra, [Aq, Aq] — 1, allows one to define a Fock space consisted of the ground state \0)g defined 
by, 

i G |0) G = 0, (5) 

and the excited states, \ti)g, generated by successively acting A G on \0)g- 

The next step is the variational procedure with respect to fie so that the energy of the Hamiltonian J2J with 
respect to the Gaussian ground state is minimized. Some of the nonperturbative effects of the Hamiltonian are 
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amalgamated to Hq by this process, which makes the correction term, Vi, to be an order of £g (defined below) smaller 
than the Gaussian Hamiltonian. The Vi term may be written as in Ref. @] as 
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where £g 



< 1 for any value of A due to the variational gap equation, 1 — £g = — In this sense, the 
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Gaussian Hamiltonian Hq describes a truncated form of H up to 0(£q). 

The final step in the variational perturbation theory is to apply the conventional perturbation method to the 
Hamiltonian J2J) by using the complete set of states, {|n)c}- In an operator method, the so called Liouville-von 
Neumann approach, this is achieved by defining the 1 st order annihilation operator, An\ = Aq + £,gB, so that the 

Hamiltonian H can be factorized as H — —— + ^(1)^(1))) U P to the first order in £g • It was shown that 

these operators satisfy the (real) q— deformed algebra 01 with deformation parameter, q 2 = 1 + £g/2- This algebra 
leads to another Fock space determined by the 1 st order creation operator, which defines the 1 st order states 

The fact that the first order perturbative correction to the Gaussian approximation is described by the q— deformed 
oscillator indicates that one may establish a better variational perturbation theory by using the q— deformed oscillator 
as the basis for the variational method. To develop the variational perturbation method based on the q— deformed 
oscillator, we start by separating the Hamiltonian Q into the q— deformed part, H qi and the perturbation term, V[: 



H = H n 



VI, 



(7) 



where V[ = H — H q , and 

Hq = '~y {a q a\ + a\a q ) 

is the Hamiltonian of the q— deformed oscillator ^l|> with a q and a q satisfying the algebra, 

it; 



[a g , at] = 1 + ea' q a q , 



(8) 
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where we use this specific form of a— algebra, since this is satisfied with the perturbative anharmonic oscillator |6|]. 
The q— ground state is defined by 

& q \0) q = 0, (10) 

and the q— excited states can be generated by successively acting d q to the ground state (|l(Jfl . 
For calculational simplicity, we introduce dimensionless operator of order 0(h°): 



H n = 



2fl n 



x m p n + p n x r ' 
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Similar operators are defined in Ref. [l5|. The explicit expression for a q , which satisfies the algebra © and make the 
Hamiltonian JSJ to be (|33|l is obtained to 0(e 2 ) in Appendix A, 
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where u* = — ( 1 e e 2 I . After inverting the expression for a q to get H® as a function of a q and a q , and taking 
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the expectation value with respect to the q— ground state (|10[l . we obtain 



(H° 2 ) q = l-e+^e 2 , 



(H$) q = 15. 



(13) 



Note that we have introduced two variational parameters Q and e. To establish the variational perturbation theory 
based on the q— Hamiltonian, we need to express the q— Hamiltonian (jHJ) in terms of the phase space variables, x and 
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p. By demanding that the q— Hamiltonian is frictionless (no px + xp term), the q— Hamiltonian may be written as, 
up to 0(e 2 ), 
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The parameters in ||SJ) and l|14l) are determined in the appendix A and 53 = . 

45 

Note that, unlike the case of the variational perturbation theory based on the Gaussian approximation, the 
q— Hamiltonian H q is not expressed in a closed form in the phase space, but is written as a series in e. Since we 
are to compute up to the 1 st order correction to the variational result, it is enough to write H q up to 0(e 2 ) as in 
Eq. d. 

The Hamiltonian H of Eq. (JTJ is then written as, H = H q + Vr, up to 0(e 2 ), where 
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with £ = — — ^- . Taking expectation value of JJJ with respect to the q— ground state Ql()[). the energy expectation value 
2\l q 

of the anharmonic oscillator becomes 
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+ 0(e 3 ). 



(16) 



The variational minimization of the energy expectation value (|16|l . with respect to f2 g and e, leads to the gap equations, 
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(17) 



3 A/i 

where £0 = — 7- The gap equations relate f to e as, 
2oj 4 



f = ^l 1 " 4 
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By using Eq. I|17[l , the ground state energy of the anharmonic oscillator becomes 
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where £(£0) and e(£o) are determined by the gap equations (|17l) . We present the values of e, £, and Eg for several 
values of the coupling, £oj m Table I. As can be seen in this table, the maximal error of the present approximation is 
less than 0.8%. Moreover, the accuracy is 200 times and 10 times better than the Gaussian and the post (including the 
2* order perturbative corrections) Gaussian approximations , respectively, for weak coupling. The maximum value 
of the dimcnsionless expansion parameter, e ~ 0.1717, is attained as £0 ~* °°j which is distinguishably small compared 
to the Gaussian expansion parameter, £g — 1, at the same limit. This clearly shows that the present method provides 
better nonperturbative information than the Gaussian approximation and 2 nd order variational perturbation results 
based on the Gaussian approximation. 



To complete the 1 st order perturbative corrections to the variational result, we need to construct the 1 st order 
creation and annihilation operators which are correct to 0(e 2 ). To do this we need to express the Hamiltonian Q as 
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TABLE I: ground state energy of anharmonic oscillator 
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a function of q— operators, a q and a). Eq. I|18fl enables one to write the coefficients in the potential Vjf, Eq. i|15[) . as 
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9i(e) = 52(e) = 3(l + ^e-^e 2 )/+^ 
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Note that 5i(e)'s are of O(l), and thus this shows the fact that V[ is order of e smaller than H q . The Hamiltonian is 
then written as, 



H = — [a q al + <r : <i„ + 
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We now want to write this Hamiltonian as a generalized deformed oscillator: 
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H = — (a {1) aj il) +a[ 1) a (1) J + 0(e 3 ), 



(22) 



[«(!))«(!)] = 1 + eaia( 1) a( 1) + e 2 a2(a( 1) a(i)) 



where the algebra defines the deformation function, F(y) = 1 + (1 + eati)y + a 2 (ey) 2 53- Since V[ is 0(e 2 ), correction 
to the annihilation operator would be of order e 2 , and thus, am can be written as 
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For this operator to satisfy the algebra (1221) . the following relations should hold for the coefficients in Eq. I|23ll : 
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«d = j = (52 -1033)^, «2 + « 2 = (52 + 10 93 ) W , «3 = - T> 
«i + uj;=0, u = YrqCsi + 6 52 + 4553), 



(23) 
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which give, to 0(e 2 ), 
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To obtain the 1 st order correction to 0(e 2 ), we may set f2 g /f2 — > 1 in Eq. (|24|) . since all the coefficients of fl q /il have 
factors of order 0(e 2 ). The 1 st order ground state is then defined by 

o (1) |0) (1) =0, (26) 

and its energy is the same as that in Eq. I|2l)|) . since the variational approximation leading to (|2(Jfl includes the 
contribution from the potential Vj, as in the case of the Gaussian variational perturbation. The n th excited states 
are given by 

|n) (1) = %4fl°>(i)> ( 2? ) 



where [n]l — [n][n — 1] • • • [1], and [n] is defined by the recurrence relation [n + 1] = with [0] = 0. The energy 

of the n th eigenstate is given by 

E n =™([n] + [n+l]). (28) 

In contrast to the ground state energy that shows a slight improvement for large £o, the energies of the excited states 
receive considerable improvements compared to those in the Gaussian and the post Gaussian approximation, since 
the expansion parameter e is smaller than that of the Gaussian approximation. 

The higher order corrections can be obtained similarly as in the above procedure. To obtain the 2 nd order correction, 
for example, we need to write the Hamiltonian (|21|l up to 0(e 3 ), and to express it as a function of and o/ny 

We then write the Hamiltonian in a factorized form, H^) — _ ^ _ (^(2)S-(2) + ^(2)^(2)) + *-^( e4 )' while a( 2 ) satisfies a 
generalized deformed algebra. 

The fact that the first order perturbative result can be expressed as a generalized deformed oscillator, Eq. (|22|) . also 
provides us with a possibility of establishing a new variational perturbation theory based on the generalized deformed 
oscillator (|2*2l . 

One of the main advantage for using the algebraic approach is that it enables one to define the thermal state and 
coherent state easily. For example, the coherent state can be written as the state |a), 



a(t)\a) = \a)a. (29) 

The construction of the state |a) from the number states can be achieved from the usual method of the q— oscillator (lif . 

It would be interesting to apply the present method to general quantum mechanical systems. For example, the 
realization of su 9 (2) by q— boson can be used as a basis of a perturbation theory for systems with spherical symmetric 
potential. More interesting would be to generalize the present method to quantum field theory by expanding the 
quantum fields in Fourier modes where each mode acts as a generalized deformed oscillator. 



This work was supported in part by Korea Research Foundation under Project number KRF-2003-005-C00010 
(H.-C.K. and J.H.Y.). 



6 



CALCULATION OF THE ANNIHILATION OPERATOR 



To have an explicit expression for the annihilation operator d qi we need the q— annihilation operator as a function 
of x and p, 



21+1 



E ^<l H 'n 



1=0 n=Q 

which satisfies Eqs. (jSJ and @. In obtaining this expression, the operator product expansion formula, 



m-\-m' — k 



k=0 



(30) 



(31) 



(32) 



is useful. The expansion coefficients, hk are given as 
ho = 1, hi = mn — m'n, 

I12 = — [m{m — 1) + m (m — 1)] nn — mm! [n(n — 1) + n (n — 1)] — 3mm' nn' , 
h 3 (m, n; m', n') = 2{f£ -/£,) + (/£ +m , - /£+"') - im'm{m - 1) ( " + ) + 3n'n(n - 1) ( m + m ' 
—3 [m(m — l)n (n — 1) — m (m — l)n(n — 1)] (m + m' — 2)(n + n' — 2), 



where /£ = 6 , 3 i . 3 

We assume that the Hamiltonian for the g-oscillator is written as the sum of the kinetic energy and potential term 
dependent only on x, 



P 



G 



(33) 



If we demand the operator a 9 to satisfy the algebra JjJJl and the Hamiltonian ||SJ| to be l|33|l . we get a q in Eq. I|12|l and 
the q— Hamiltonian 



H a = 



The parameters in JHJ) and 1)14(1 are determined as 

n = 



1 1 ,\ 4efi 3 

23 
~45' 
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36 2 53 ^ 
= ro— -, 53 



(34) 



(35) 



/t 2 



The present result can be easily checked using mathematica, in which the noncommutative multiplication is imple- 
mented. 
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